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Dilepton production from the quark-gluon plasma (QGP) phase of 
ultra-relativistic heavy-ion collisions is computed using the leading-order 
(3+l)-dimensional anisotropic hydrodynamics. It is shown that high-energy 
dilepton spectrum is sensitive to the initial local-rest-frame momentum- 
space anisotropy of the QGP. Our findings suggest that it may be possi¬ 
ble to constrain the early-time momentum-space anisotropy in relativistic 
heavy-ion collisions using high-energy dilepton yields. 
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1. Introduction 

The enormous amount of data collected in ultra-relativistic heavy-ion 
collision experiments at the RHIC (Relativistic Heavy-Ion Collider) and 
LHC (Large Hadron Collider) have indicated that a new state of nuclear 
matter called quark-gluon plasma (QGP) is created in these events. The 
data indicated that the QGP behaves like an almost perfect fluid, which 
may be, to a great extent, described within the framework of relativistic 
dissipative fluid dynamics (see the recent review [1]). Recently, the study 
of the mechanisms leading to, and the level of momentum-space isotropiza- 
tion/thermalization of partons comprising the QGP became an active re¬ 
search area [2, 3]. Understanding of these phenomena is crucial not only for 
setting up the initial conditions for fluid dynamical frameworks, but also for 
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judging applicability of the fluid dynamical approaches to the early-time dy¬ 
namics of high-energy heavy-ion collisions. Various theoretical approaches 
have attempted to explain the mechanisms driving the QGP towards the lo¬ 
cal equilibrium state in both the strong and weak coupling limits, including 
gauge-gravity (AdS/CFT) duality, perturbative quantum chromodynamics 
(pQCD), and saturated initial-state color glass condensate (CGC) models 
[3]. These results suggest that during the early stages of the evolution the 
matter possesses substantial local momentum-space anisotropies, with the 
transverse pressure significantly exceeding the longitudinal pressure. In or¬ 
der to test these findings, one would like to find experimental observables 
that are sensitive to the early-time degree of momentum-space (an-)isotropy 
of the QGP. 

High-energy (E > 2 GeV) electromagnetic probes, in particular dilep¬ 
tons (produced from the decay of virtual fotons) and real photons, seem to 
be perfect observables for this purpose, since they are produced mainly in 
the early stages/central region of the collision and they are weakly coupled 
to the QGP. This means that they leave the system almost undisturbed (for 
recent reviews see [4, 5]). In order to study the effects of local-rest-frame 
(LRF) anisotropies in the system, the calculation of these observables has 
to be convoluted with the fluid dynamical framework which treats these ef¬ 
fects in a reliable manner. Along these lines, in this work, we review recent 
results obtained by folding the dilepton production from the momentum- 
space anisotropic QGP with the QGP background evolution obtained us¬ 
ing the recently developed framework of (3+l)D leading-order spheroidal 
anisotropic hydrodynamics [6, 7, 8, 9, 10, 11, 12, 13, 14], As reported orig¬ 
inally in Ref. [15], we find that the high-energy spectrum of dileptons is 
quite sensitive to the early-time momentum-space anisotropy of the QGP, 
and, therefore, this observable may potentially be used to experimentally 
constrain the degree of early-time momentum-space anisotropy. 


2. Dilepton rate in an anisotropic plasma 

At leading order in the electromagnetic coupling, 0(a 2 ), the dilepton 
emission rate comes from quark-antiquark annihilation. The calculation 
of the rate in an anisotropic QGP was performed using relativistic kinetic 
theory in Refs. [16, 17]. Therein, it was shown that the production rate is 
given by 


dR l+l 

d 4 P 


^lj(cOs0 Pl )j a - f q (pi y4 +£ COS 2 0 pi , A^ 
xfq(^\J{E-pi) 2 + ^{piCOS0 pi - Pcosdp) 2 , a) , 


(1) 
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where 


a± 


M 2 

2(E ~ Pcos(6 P ± 6 Pl )) ’ 


(2) 


and 


X = (2piP sin 0psin 0 Pl ) 2 — [2pi(E — P cosOp cosd Pl ) — M 2 } 2 . (3) 

Above it was assumed that the three-momenta of the quark and the lepton 
pair are parametrized using spherical coordinates in the following way 

Pi = Pi(sin0 pi cos </> pi , sin 0 pi sin <^ pi , cos 0 pi ), 

P = P(sin$pcos(?ip,sin#psin(/>p,cos0p), (4) 

where the z-axis is oriented along the direction of anisotropy n (which is 
taken to point along the beam line direction), and E and M are the energy 
and the invariant mass of the lepton pair, respectively. The phase-space dis¬ 
tribution functions of the quarks and antiquarks, which have to be specified 
in (1) are assumed to be equal and, at leading order, given by the spheroidal 
Romatschke-Strickland (RS) ansatz [18] 

fq{q) (p, A ) = fq°q)(Vp 2 + £(P ' n) 2 , A) . (5) 

In Eq. (5) the anisotropy parameter £ describes the strength and type of 
anisotropy. The physical situations suggested by microscopic models dis¬ 
cussed in Sec. 1 correspond to £ > 0 (oblate distribution), however, we allow 
£ to be anywhere in the range — 1 < £ < oo. The parameter A is the trans¬ 
verse momentum scale, which coincides with temperature in the case of local 
isotropic equilibrium (£ = 0). The isotropic distribution is assumed to 
be a Fermi-Dirac distribution. In general, the transverse momentum scale 
and the anisotropy parameter may be arbitrary functions of space-time 
point, which in Milne parametrization reads AT = (rcosh^, x j_, t sinh <f) . 
Herein, we use standard notation for longitudinal proper time r = Vt 2 — z 2 
and the space-time rapidity ? = tanh ~ 1 (z/t). In this study, the functions 
£(X) and A(X) are taken from from fluid dynamical simulations performed 
using the anisotropic hydrodynamics framework, see Sec. 4. 


3. Dilepton spectra 

With the dilepton emission rate in hand, the final invariant mass and 
transverse momentum spectra may be obtained by performing the following 
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phase-space integrations 


dN l+l 

MdMdy 

dN l+l ~ 

P±dp±dy 



M max 2ir /■ + , s 

y Mmin MdM J o d( t>v j dA x—^r(i{x), A(xy), 


(6a) 

(6b) 


respectively. The four-momentum of the dilepton pair is parametrized in 
the following way 


P^ = (mj_ cosh y, p± cos 4> p ,p± sin 4> p , m± sinh y) . (7) 

Since we work in the Milne coordinates the space-time measure is d^X = 
rdr dsd 2 x±. Before evaluating formulas (6a) and (6b) one has to remember 
to transform the LAB momentum of the dilepton pair to LRF (where Eq. (1) 
was obtained) using p'^ = h!‘ v p v . In this case, the general Lorentz boost 
tensor depends on the four-velocity of the local rest frame h! L v K(A)). 


4. Hydrodynamic evolution 

Although the system studied here is strongly anisotropic, it is assumed 
that it behaves collectively, and may be effectively described using a finite 
set of macroscopic degrees of freedom using a fluid dynamical framework. 
For this purpose, we use the anisotropic hydrodynamics approach. 

At leading order, the evolution equations of anisotropic hydrodynam¬ 
ics may be derived by taking moments of the relativistic Boltzmann ki¬ 
netic equation in the relaxation time approximation and using the RS 
ansatz (5) for the distribution function. In this way, the resulting energy- 
momentum conservation equations = 0 and particle production equa¬ 
tion = u p (A r eq — N ^) /r eq lead to five partial differential equations for 

five functions of space-time, i.e, anisotropy parameter £, transverse momen¬ 
tum scale A, and three independent components of four-velocity u fl . The 
relaxation time is given through the relation r eq = 5fj/(2T), where fj = rj/s, 
with T] and s being the shear viscosity and entropy density, respectively. The 
initial condition for hydrodynamic evolution of A is specified on a constant 
proper time surface using a mixed optical Glauber model. We also assume 
that the initial flow in the transverse direction is vanishing, while in the 
longitudinal direction one initially has Bjorken flow, v z = z/t. For details 
of the hydrodynamical calculations and values of the parameters used, we 
refer the reader to Sec. 4 of Ref. [15]. 






dN/MdMdy [GeV 
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5. Results 

The hydrodynamic model presented in previous Section is used to sim¬ 
ulate a (3+l)-dimensional minimum-bias lead-lead collision at LHC beam 
energy of 2.76 GeV. After the hydrodynamic evolution is determined the 
dilepton spectra are computed using formulas (6a) and (6b), assuming that 
the emission stops when the temperature drops below a critical temperature 
of 175 MeV. In the calculations, we restrict ourselves to the high-energy part 
of the spectrum by setting the cuts to = 1 GeV and p™ ax = 20 GeV 
for transverse momentum integration, and M mm = 1 GeV and M max = 20 
GeV for the invariant mass integration. In Fig. 1, we present exemplar re¬ 
sults for the invariant mass (left) and transverse momentum (right) spectra 
for various values of the initial momentum-space anisotropy. One observes 
that the transverse momentum spectra in particular is extremely sensitive 
to the initial anisotropy assumed. We observe that, with increasing ini¬ 
tial anisotropy, the spectra becomes flatter. We checked that the observed 
effects are larger for 77 /s > 1/47r. 




M [GeV] [GeV] 

Fig. 1. (Color online) The minimum-bias invariant mass (left) and transverse mo¬ 
mentum (right) spectra of dilepton pairs calculated within anisotropic hydrody¬ 
namics for various values of the initial anisotropy £0 and the AdS/CFT lower 
bound value of shear viscosity to entropy density ratio, p/s = l/(47r). 


6. Conclusions 

Based on the model results obtained, we find that the invariant mass 
and transverse momentum spectra are sensitive to the initial anisotropy of 
the QGP. Our findings provide possibility to constrain the level of isotropi- 
sation of the strongly interacting medium created in relativistic heavy-ion 
collisions by measuring the high-energy part of the dilepton spectra. It 
would be very interesting to confront our findings with the experimental 
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data including additional dilepton emission from the hadronic phase. More¬ 
over, it would be interesting to go beyond the leading-order formulation of 
the anisotropic hydrodynamics by including NLO corrections [19, 20], or to 
use more complete leading-order formulations of anisotropic hydrodynamics 
[ 21 , 22 ], 
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